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Abstract 

We present and solve a continuum theory of vacancy-mediated diffusion (as 
evidenced, for example, in the vacancy driven motion of tracers in crystals). 
Results are obtained for all spatial dimensions, and reveal the strongly non- 
gaussian nature of the tracer fluctuations. In integer dimensions, our results 
are in complete agreement with those from previous exact lattice calculations. 
We also extend our model to describe the vacancy-driven fluctuations of a 
slaved flux line. 
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I. INTRODUCTION 



Random walks and associated diffusion processes are ubiquitous in Nature 0. Our 
understanding of their properties has grown steadily over the years, with contributions from 
workers in many fields, but there are always surprises in store as we learn to ask new and 
more probing questions. 

An interesting example of a non-trivial random walk problem is that of vacancy-mediated 
diffusion (VMD) This applies to the case of a crystal containing a low density of 

vacancies; the question being, what are the fluctuations of a tagged particle (or impurity) 
which moves only due to exchange with wandering vacancies? This process has a more 
generic application as it is one of the simplest 'slaved diffusion processes'. It has been 
studied in various guises over the years || and a lattice formulation was solved exactly 
by Brummelhuis and Hilhorst |4[] in the late Eighties. The most revealing aspect of their 
solution is the strongly non-gaussian nature of the tagged particle's fluctuations. The same 
lattice model received alternative exact treatments recently ||. 

Given one is generally interested in the scaling behaviour of such systems, it is useful 
to have a more coarse-grained treatment to hand which allows one to access the (hopefully 
universal) long-wavelength behaviour directly, without having to resort to a sophisticated 
lattice calculation. A well-known example of such a situation is the critical behaviour of a 
ferromagnet || , which may be studied within the context of the lattice (Ising) model (where 
an exact treatment is only possible in dimensions d = 1, 2), or at a more coarse-grained level 
via the 4 field theory (where controlled calculations can estimate the exponents in the most 
interesting case d — 3) 0. 

The purpose of this paper is to construct and analyse a continuum (or coarse-grained) 
theory of vacancy-mediated diffusion. Fortunately the continuum theory is exactly solvable 
in all dimensions, which allows a direct comparison with the previous lattice calculations 
Q - complete agreement is found. Given this situation, one has confidence in applying the 
continuum theory to situations where a lattice formulation would be more difficult, if not 
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intractable. We consider one such situation here; namely the motion of a directed (flux) line 
under the action of diffusing vacancies. 

The outline of the paper is as follows. In the next section we define the lattice model of 
vacancy-mediated diffusion more precisely, and, using this as a base, construct the continuum 
theory. We then highlight the 'mesoscopic' physics contained within this coarse-grained 
model. In Section III we analyse the continuum model at the level of mean-field theory 
(MFT). The apparent physics within MFT is simple, and the solutions of the mean- field 
equation yield results which contain some correct scaling information (i.e. the length-time 
scaling), but completely miss the more interesting statistical aspects; namely the strongly 
non-gaussian nature of the fluctuations as found from the lattice calculations. Thus we 
are led to attempt an exact solution of the continuum theory. This is made possible by 
formulating an infinite-order perturbation theory, as described in Section IV. The results 
for various dimensions are derived and presented in Section V - the non-gaussian statistics 
are seen to be completely reproduced, with the added advantage that the results from the 
continuum theory are valid for arbitrary dimension d. In Section VI we give an illustration 
of the utility of this coarse-grained approach. We formulate a continuum model for the 
transverse fluctuations of a d + 1 directed (flux) line due to VMD of pinning centres. We 
end the paper with Section VII which contains our conclusions, along with some suggested 
extensions of the present work. 

II. FORMULATION OF THE CONTINUUM THEORY 

We first consider a simple lattice formulation of VMD with which we shall motivate our 
continuum theory. To be specific, let us consider the following lattice model. Each site 
r of a rf-dimensional hypercubic lattice contains a spin S r which may take the values ±1. 
[In the crystal lattice application one would take all spins as 'up', thus referring to crystal 
atoms, except for the tagged particle (or impurity) which would be assigned a 'down' spin.] 
A spin may only alter its value when involved in an exchange with a single diffusing vacancy 
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located at the position R(t). As only the position of the vacancy is relevant we may give 
the vacancy a spin +1 for notational convenience. [We shall assume for simplicity that only 
one vacancy is present - the generalization to many vacancies is straightforward.] There are 
various ways to describe this model. For instance, one can define a probability distribution 
for the vacancy position and the position of the 'up' and 'down' spins, and then write a 
master equation for its evolution. Alternatively, one can regard the process as a stochastic 
cellular automata (SCA), and attempt to write explicit SCA rules for its operation. We 
shall follow the latter approach. 

The rule for updating the vacancy position is written as 

R(t + St) = R(t) + l(t) , (1) 

where l(t) is a unit lattice vector drawn with equal probability from the 2c? possible choices. 
The update of the spins is written down as follows: a spin will remain unchanged unless it 
is located either at the vacancy position at time t, or at the subsequent vacancy position at 
time t + St. Thus we have 

S r (t + St) = S r (t) + 5 rj R( t ) [S-R( t +st) (t) - S r (t)] + Gv,R.(t+«) [Sn(t) (t) - S T (t)] , (2) 

which may be re-arranged (with the help of ([!])) in the appealing form 

S r (t + St) — S r (t) = [5 , R (j) + i(t)(t) — SR(t)(£)] [S r ,n(t) — £r,R(t)+i(t)] • (3) 

The factorization of the SCA update for S T (t) gives us reason to hope that a simple 
Taylor expansion may be invoked to yield an accurate continuum limit, since the leading 
order terms will appear in a multiplicative fashion (as opposed to terms occurring in an 
additive fashion, in which case one is unsure of their relative importance). Thus we shall 
take the simplest continuum limit: First, the position of the vacancy is described by a real 
vector R(t) 6 1Z satisfying 

dR/dt = f (t) , (4) 
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which is the continuum equivalent of Eq. ([!]). The random lattice vector 1 has given way 
to a random vector £ (with an implicit factor of y5t) which is drawn from a gaussian 
distribution P[£] with zero mean, and covariance = D5 a ^5(t — t'). [Henceforth, 

angle brackets denote an average over P.] Second, the spin variable S r (t) will be replaced 
by a coarse-grained 'magnetization density', or 'order parameter' (OP) denoted by the field 
4>(r,t). Expanding each term of Eq.(||D to first order, we obtain 

d(f>(R(t),ty 



dt<f>(r,t) = -A 



d t 6 d (T-TL(t)) , (5) 



dR(t) 

where A is a phenomenological parameter with dimensions L d T. This last equation represents 
our continuum theory of VMD. It loosely resembles a Langevin equation, but the extremely 
implicit appearances of the noise £ throughout the equation forbid such a simple designation. 
For instance, it is not clear how one would write down a dynamical equation for the OP 
probability distribution P[(f>, t] - it would certainly not fit within the standard Fokker-Planck 
category. However, taxonomy is of little importance to us. We shall accept the equation 
on its own merits. As a first step in this direction, let us probe its more obvious physical 
content. 

We can ascribe an independent meaning to each of the two factors on the right-hand-side 
of Eq.([|). The second of the two is simple enough in spirit - it allows temporal change of 
the OP only in the neighbourhood of the vacancy position R(t). The first factor guarantees 
(via the directional derivative) that the OP temporally changes only when the vacancy 
moves through a region in which the OP has spatial variation. Furthermore, the amount of 
temporal variation is linearly coupled to the amount of spatial variation (with a strength A). 
This appears reasonable with regard to the physical properties of VMD we wish to model. 
Whether one would have written down such a continuum model a priori based on these 
considerations is unclear. However, given the simple 'derivation' from the SCA, one finds 
the a posteriori physical motivations satisfactory. 

It is useful to note that analytic progress on Eq.(^) is very difficult without first recasting 
it in Fourier space (so as to rid ourselves of the implicit nature of the noise). Denoting the 



Fourier transform (FT) of the OP by <f)(k, t) we have 

d t $(k,t) = -X G(k,t) J dki G(k u t)*4>(k u t) , (6) 

where 

(5(k,t) = (£-k)exp[ik-R(t)] , (7) 

and dk = d d k/(2n) d . 

We shall now discuss the initial conditions. As for the vacancy, we simply need to define 
its initial position Ro = R(t = 0). The initial condition 0o( r ) for the OP will vary according 
to the physical system we are interested in modelling. As regards the case of a tagged 
particle (or impurity) being driven by the wandering vacancy, we take the initial condition 
on the lattice to be all spins up, bar one down spin (representing the tag) located at some 
point (the origin, say). In the continuum, this may be described by 

0o(r) — A- B8 d {r) . (8) 

Alternative scenarios may be investigated by modifying the initial OP distribution. For 
instance, setting <J)q to be a step function would be appropriate for modelling the vacancy 
mediated roughening of an initially straight domain wall [H. 

Once the initial conditions are defined we may formally integrate the equations of motion 
to give 

t 

R(t) = R + J dt' f (*') (9) 
o 

and 

t 

4>(k, t) = o (k) -\Jdt' G(k, f) J dk x G(ki, t'Y 4>{ki, f) . (10) 
o 

This completes our formulation of the model. The main focus of this paper is to calculate the 
mean OP density p(r, t) = (4>(r,t)), which is the quantity analogous to the tagged particle 
distribution function calculated in previous lattice studies. In the next section, we shall 
present a brief mean-field analysis of p, whilst sections IV and V contain a description of 
the exact solution for p from Eqs.(f|) and (|10|) above. 
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III. MEAN FIELD THEORY 



The purpose of this section is to indicate how much one can learn about the system from 
a simple, yet uncontrolled, mean field theory (MFT). We shall find that MFT predicts the 
correct length-time scaling, but misses the non-gaussian nature of the fluctuations in VMD. 
This error persists in all dimensions (bar d — 0). 

We shall define the MFT to be used here in an operational sense. Namely, we perform 
the simplest possible average of the OP equation of motion (P), by replacing the average of 
the right-hand-side by the product of two separate averages. Explicitly, we have 

d t p(k,t) = -X J dh (GfaQGfat)*) p(k 1; t) . (11) 

We refer the reader to appendix A for the evaluation of (G(k, t)G(ki, £)*}. The averaging 
necessarily introduces a temporal cut-off to (which sets the implicit correlation scale of the 
white noise £). We shall only ever work to leading order in 1/to- Eq.flTTD now takes the form 

d t p(k,t) = ———k a [ dh fc« p(ki,t) e-^m^-^t+i^yno ? (12) 

where the momentum components are indicated by a superscript for future notational con- 
venience. 

This integral equation may be recast in a more illuminating form by inverse Fourier 
transforming the density. After some re-arranging, one has 

d tP (r,t) = ~- J dh fc« e- ikl - r p(k!,t) J dk (k a + ^) e -P/2)^-*k.(r-R ) _ ( 13) 

Now the inner integral may be written as (id r a + k")g(r — H ,t). The Green function 
g(r, t) = (2nDt)~ d / 2 exp[— r 2 /2Dt] is the probability density for the vacancy (in other words, 
g is the solution of the Fokker-Planck equation corresponding to Eq.(^)). One may then 
manipulate the outer integral over ki in terms of the inverse FT of the density to obtain the 
following partial differential equation for p(r, t) : 

d tP =^-V-[g(r-R ,t)V] p . (14) 
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As a MFT, the above equation makes good sense. It contains the simple physical in- 
formation that the OP density undergoes a diffusion process, but with the twist that the 
diffusivity is not constant, but proportional to the probability density of the wandering 
vacancy. One would have had little trouble in writing down such a MFT using physical 
arguments alone. 

A complete analytic treatment of Eq.fll4|) is beyond our present remit. However, the 
limit r 2 <C Dt is trivially solved, since to leading order the equation reduces to 

e(t)d t p = ^V 2 p , (15) 

where e(t) = {t Q /2\)(2nDt) d / 2 . One may solve the above equation using FT, and with the 

initial condition specified in Eq.® we obtain p = A — Bg(r, t(£)), where r(£) is the effective 

t 

time scale of the tagged particle, and is given by r(t) = J dt' e(t') -1 which, ignoring 

to 

numerical prefactors, takes the form 



r(t) ~ < 



TtMt/to) d = 2 (16) 



A 



d > 2 



(Dt ) d / 2 

This solution firstly tells us that the OP density always has a gaussian envelope (at least 
for r 2 <C Dt, which encompasses most of the physically interesting scales, since t ^> r), but 
that the temporal spreading of the envelope is not that of a random walker (Ar ~ 
but much reduced. In two dimensions the spreading increases only logarithmically in time 
(as the vacancy's random walk is only just recurrent), and for d > 2 the spreading halts 
altogether after some finite time, which is accounted for by the vacancy having 'fled the 
scene', never to return. 

The main interest in VMD is not so much in the renormalized time scale, which one 
can argue for on simple physical grounds, but rather in the non-gaussian nature of the OP 
fluctuations. These non-gaussian fluctuations were exactly calculated in lattice theories, and 
were generally found to have tails which decay slower than a gaussian. It is the aim of the 
following two sections to reproduce these features from the continuum theory. 
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IV. SOLUTION VIA PERTURBATION THEORY 



We have seen the relative failure of MFT, which is essentially due to imposing a strict 
locality on the OP equation of motion. In fact, the evolution of the OP is non-local as 
seen from Eq. (|5]) . A more systematic treatment is required to handle the subtle correlations 
between vacancy and OP. Given the functional nature of Eq. ([)]), the most useful analytic 
technique would appear to be perturbation theory. Since the equation of motion is linear, 
we shall not encounter an exponentially divergent number of terms at higher orders; rather, 
each order will contain only a single term. 

So referring to the time integrated (and Fourier transformed) evolution equation for the 
OP, namely Eq.([lfj), we make the substitution 



4>(k,t) = J2^Xn(Kt) , (17) 

n=0 

where Xo(k, t) = 0o(k). Equating powers of the coupling A yields (for n > 0) 

t 

Xn(k, t) = -J dif G(k, t') J dk' G(k', 0*Xn-i(k', t') . (18) 



This relation may be iterated to give the explicit solution for each order of perturbation 
theory as 

t tn-l r- n _ 1 

X„(M) = (-1)" /<*<!••• J dtnG^tt) J] J dk m G(k m ,t m yG(k rn ,t m+1 ) 
Lm=l 

x J dk n G(k n ,t n yMK) . (19) 

In principle, the solution given above may be used to calculate a range of spatio-temporal 
OP correlation functions. Our present aim is more modest - we shall perform a direct average 
of each term in the perturbation series in order to obtain the OP density p(r,t). We have 

oo 

p(M)=0o(k) + £^<Xn( M )> (2Q) 
n=l 

Using the explicit form for the G functions as given in Eq.(|7|) we may write 
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t 




"n— 1 ■ 


J dtf 


' j di n k ^ 


n/ 








.m=l 



Ji, h,Pmh.a m 



(21) 

where the function Q n represents the following average: 

Q„({k m ,t m ;a m ,/3 m }) = (-l) n ( ft r m ^ m exp \iR(t m ) ■ (k m _ x - k m )] \ , 



(22) 



\m=l 



with the time ordering t\ > t 2 > ■ ■ ■ > t„, and the notation k = k. At this stage of the 
calculation we can see clearly the remaining steps. First, we must perform a multivariate 
average in order to determine the function Q n . Second, we must perform the n-fold integral 
over the momenta {k m }. Third, we must perform the n-fold integral over the intermediate 
times {t m }. Finally, we are left to resum the functions \ n which will yield the FT of the mean 
OP density. In the absence of hindsight, it is somewhat surprising that all these steps may 
be performed exactly for arbitrary dimension d. The remainder of this section will consist 
of the details of the first two steps, whilst the third and fourth steps will be presented in 
the subsequent section as they are dimension-specific. Henceforth, we shall take the initial 
position of the vacancy to be at the origin: Ro = 0. This slight loss of generality (which is of 
no physical significance in the long-time regime) is more than compensated by calculational 
simplicity. 

We begin with the first step; that of determining Q n . We refer the reader to Appendix 
B in which this multivariate average is evaluated. The result is 



Qn({k m , t m , Ol m , 

/U) 



tr 



n [* 



- Dt (k am - k^)(k 0m - k^) 



m=l 



x exp 



Tr(k k m ) (t m t 



m+l / 



(23) 



where the symbol t n+1 = 0. 

To proceed with the momentum integrals, it is necessary to decide upon an initial condi- 
tion. We shall use that given in Eq.(|J), since our interest is presently focussed on the vacancy- 
driven diffusion of a tagged particle. The generic momentum integral (for 1 < m < n — 1) 
takes the form: 
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J dk m kkkzr* [5 am , Pm - Dt (k a - - ex**- - fcft-) 



exp 



2 (k km) (^m ^rri+l / 



= [D(t m - t m+l )}- l [2nD(t m - Wi)r J [^ w + ^(^ - Wi)^ m ^ m+1 J , (24) 

where we have worked to leading order in (1/io)- The final momentum integral involves the 
OP initial condition which in Fourier space takes the form 0o(k) = {2n) d A5 d — B. Its 
evaluation yields 



dk n ktMK)\s anA - Dt (k a « - K n ){k^ - kt) 



exp 



^"(k k n ) t n 



-Bk an (2tt Dt n y d/2 , (25) 



to leading order in (l/t Q ). Collecting our results from Eqs.(£3|), fl24|) and (|25|) , and substi- 
tuting back into (^l|), we have 



t r tl tn—1 

(x n (k,t)) = -BD[-(27iDyt ]- n Jdt 1 J dt 2 {t x -t 2 )-^--- J dt n (t n _ X -,„> 

^0 

xk ai \l[[S* m ,a m+1 + D(t m - t m+1 )k am k a ^]\k a ™ , (26) 

I m=l ) 

where we have introduced 7 = d/2. An important point must be mentioned at this stage. 
The intermediate time integrals above appear to be divergent at one or both of their lower 
and upper limits. This divergence is regularized by the existence of the microscopic time 
scale to- As mentioned before, this scale appears as an effective correlation time in the white 
noise process. It may be taken to be arbitrarily small (with respect to any 'experimental' 
time scale in which one is interested). Therefore, any time integral limit is naturally softened 
by to. 

We have now completed two of the four steps in the calculation of p(k, t). To proceed 
further requires the evaluation of the n-fold integral over the intermediate times which is 
sensitive to spatial dimension, and is presented in the next section. 
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V. MEAN OP DENSITY IN VARIOUS DIMENSIONS 



As mentioned above, we must now specify the spatial dimension of interest. In fact there 
are three cases: i) d > 2, ii) d — 2, and iii) d < 2. These cases were already apparent 
within MFT, and arise because of the qualitatively different behaviour of the vacancy's 
random walk. In case i), the vacancy will essentially 'disappear' from the vicinity of the 
tagged particle after a finite time. In case ii), the vacancy's walk is marginally recurrent 
and we expect a slow, but steady, evolution of the OP density for all times. In case iii), the 
vacancy's walk is 'strongly' recurrent, and thus the cross-section of vacancy-tag collisions is 
always 'large'. It is the aim of this final stage of the calculation to replace these qualitative 
descriptions by precise results. We shall analyse the three cases in turn. 

A. d > 2 

In this and the following two subsections we shall take advantage of the smallness of 
to- Integrals which are apparently divergent will be regularized using t , and only the most 
singular contribution will be retained. We stress that this form of regularization is not a 
mathematical manoeuvre, but is entirely consistent with the physical meaning of white noise; 
namely a noise process which is a limiting form of a microscopic process with a correlation 
time t . 

Referring to Eq. fl2~6f) we adopt the following strategy to evaluate the integrals. First, 
we explicitly contract the n-fold momentum product, which will yield 2 n_1 terms in n sets; 
terms in the m th set being characterized by a factor of fc 2 ( m+1 ) (where m counts from to 
n — 1). A term in the m th set will also carry a string composed of m different factors of 
the form (tj — t J+ i). For each term, we perform the time integrals in order, starting from 
t n , keeping only the most singular contribution at each step, and being careful to include 
the appropriate time-difference factors in the numerator (from the string). This procedure 
is fairly simple for d > 2, since the only integrals one encounters are 
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One finds that each term of a given set has the same value after integration. In detail, 
the m th set contains C^ -1 equal terms of value 

/ -. \ n—m—l / ~ \ m / ^ \ 



, 7 -l 



V(7-WV V(7"W 
Thus, the series composed of the n sets is nothing more than a binomial series; which is 
trivially summed. The dominant contribution from the n-fold time integral may therefore 
be combined with the constant prefactor of Eq. (|26|) to give 



(X»(M)> 



BDk 2 



(7-1X27^0)7 



(7 



^Dk% 



1 ^ n-1 



(28) 



This ends the third step (namely the intermediate time integrals). The last step is to 
reconstruct p from the infinite sum given in Eq. (p0|) , using Eq.(^) above. In the present 
case, this sum is seen to be simple, as the series (in powers of A) is geometric. Explicitly 
evaluating the sum yields the final result in the form 



p(k,t) = A(27r) d 5 d (k 

where the momentum scale is given by 

(7-1) 



B 



1 + {k/Kf 



(29) 



A* 



2Dtr 



A 



(30) 



It is of interest to recast this result in real space. For instance, in the physically pertinent 
case of d = 3 one has 



p( r , t) — A — 



B 



5\v) + 



4nr 



-Ar 



(31) 
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Thus, for d > 2 we find that the initial 5-function of the OP density is evolved such that 
only half of its weight is smeared. That half which is smeared attains a Lorentzian profile 
in Fourier space, with a momentum scale A as given in Eq. (|30|) . This scale is seen to be 
an effective UV cutoff (1/y/Dto) renormalized by the vacancy-tag coupling A. Note that 
this scale is independent of time, in accordance with our expectations. Note also that the 
smearing of the OP density is strongly non-gaussian. In real space, the smearing creates an 
OP density which is exponential in form, as seen explicitly for d = 3 above. 

One final point: we expect that as d increases, less of the initial 5-function (in the OP 
density) will be smeared, since the vacancy will disappear from its vicinity with increasing 
efficiency The fact that we find exactly half of the 5-function to be smeared, for all d, is a 
consequence of starting the vacancy's walk precisely at the location of the (^-function. Had 
we chosen R 7^ 0, the <i-dependence of the 'smearing fraction' would have been apparent. 



B. d = 2 

We now turn to the marginal case of d — 2. Within MFT we found that the root-mean- 
square fluctuations of the tag (i.e. the smearing of the initial 5-function in the OP density) 
grow as [ln(t)] 1//2 . We expect this slow growth to be retained within the exact solution. Our 
main interest is in how the functional form of the OP density differs from the gaussian found 
in MFT. 

In exactly two dimensions, the n th -order contribution to the OP density, as given in 
Eq.(f26|), takes the form 





x/ 



r 1 



(Xn(k, t)) = -BD(-2nDt )- n J dh J dt 2 (h - t 2 )~ 2 ■■■ j dt n {t n ^ - t n ) 

L 

W nVcw*^ +D(t m -t m+l )k a ™k a ^])k a « . (32) 

lm=l J 



Our strategy for evaluating the time integrals is the same as before. We multiply out the 
integrand to form a total of 2 n ~ 1 terms arranged in n sets. Each term is integrated over the 
n intermediate times, with only the most singular piece retained from each integral. Care is 
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taken to include the appropriate factors, for a given term, in the numerator. The integrals 
one encounters are given below (with p > 0): 

t-to 



/ 

to 

t-t 



(t-s) 2 s 



tot 



ds 



to 

t-tQ 



ds 



[Hs/t )r _ [Ht/ to )r 

> 

Kg/to)g _ (p + 2) [ln(t/t )]P +1 
(t-s)s ~ (p+1) t 

MV*o)P 1 



(33) 



to 



As before, one finds that each term of a given set has the same value after integration. 
In this case the m th set contains C^ _1 equal terms of value 

D m k 2{m+X) ( x /f )n-m-l ^ /t Q )]™+l . 

The n sets form a binomial series which is trivially summed. We find 



(x»(M)> 



5DA; 2 ln(t/t 



1 



n-1 



2?T.Dt§ 



1+DkHo \n(t/t c 



(34) 



27rDt 

The final step is to sum over the functions (% n ) as prescribed by Eq. (pU|) . As before, this 
series is geometric and the sum may be immediately performed to give 

1 



p(k,t) = A(27r) d 5 d (k) - B 



where the momentum scale is given by 



l + (A;/A) 2 ln(t/t c 



(35) 



Dt 



1 + 



A 



(36) 



Inverting the FT yields our final result 



p(r,t)=A 



BA 2 



Kn 



Ar 



[ln(t/t )]V2 



(37) 



2vr \n(t/t ) 

where K is the modified Bessel function of zeroth order ||. 

We see that the gaussian envelope for the spreading of the OP density (as found in 
MFT), has given way to a completely different form, namely the Bessel function Ko- This 
is in complete agreement with the previous exact lattice calculations Q. 
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C. d< 2 



Finally we consider VMD for d < 2. Within lattice calculations the case of d = 1 may 
be studied either within the geometry of a chain, or within the geometry of a strip of finite 
width. In the former case, the situation is trivial, since the tag can only be moved back 
and forth to one of two sites as the vacancy passes by. In the latter, the smearing of the 
tag distribution function is non-trivial; the gaussian fluctuations of MFT giving way to a 
stretched exponential. In this subsection, we shall see how to recover these results, along 
with their generalization to arbitrary d G [0, 2]. 

Our starting point is the expression for \ n given in Eq. ([26]) . We shall adopt the same 
strategy as before to evaluate the n-fold integral over the intermediate times. In this case 
we encounter the following integrals (with p > 0) 

ds 



to 

t-to 



f ds 1 „ B(l- 7 ,(p + l)(l- 7 )) 

J (t- S )7 S (P+1)7-P ~~ t (p+2) 7 -(p+l) ' 

to V ' 

1 t (p+l)(l- 7 ) 



ds 



to 



s (p+lfr-p (p + l)(l_ 7 ) ' 

where B(a,b) is the Beta function |J. 

In the two previous subsections, we were able to extract the most singular contributions 
from the 2 n_1 terms in the integrand of Eq.(^), and we found that these contributions 
formed a binomial series which was then easily summed. In the present case, this simple 
summability is lost, due to the presence of the Beta functions. However, we retain the 
feature that each of the C^j -1 terms within the m th set are equal in value. Thus, after some 
manipulations, we may reduce the function (x n ) to the form 

{Xnl ' n [-7«o(2 I rDt„)1» h ' r(l + (m+l)(l- 7 )) ' 1 ' 

In order to perform the summation over m, we introduce Hankel's representation of the 
Gamma function 0] : 
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where the contour C runs from minus infinity above the negative real axis, encircles the origin 
clockwise, and then returns to minus infinity below the real axis. Using this representation 
for the Gamma function appearing in the denominator of the summand of Eq. (|39|) , we may 
explicitly perform the (binomial) sum. Each function (x n ) now takes the form of an integral 
over r, with n appearing only as a simple power in the integrand. Thus the sum over these 
functions (as dictated by Eq.([H])) is again geometric and may be performed with ease. One 
then has the following integral expression for the FT of the mean OP density, valid for 
< d < 2 (i.e. < 7 < 1): 

1 



p(k,t) = A(2n) d 5 d (k) - B^- f dr- 
where the renormalized UV cutoff is given by 



r 1 -r + (A;/A) 2 (t/fo) 1 -T 



(41) 



A 



7 r(i- T ) J Dt c 



i + 



7(27LPt )Ttc 

A 



(42) 



We now wish to extract the scaling behaviour of the mean OP density from the above 
expression. For convenience we define 5p = A — p, which is initially a 5-function with 
amplitude B. Let us first specialize to d = 1. In this case one may simplify the above 
integral considerably, using a procedure outlined in Appendix C. The result is the following 
scaling function: 

B z r ds 



6p{T,t) = — ~ f -^exp[-s 2 -z 2 /4s] , 
7i r J Js 



(43) 

where the scaling variable is z = rA/(4t/to) 1 ^ 4 - This integral is easily analyzed for both 
z < 1 and z> 1. In the former case, one finds 

B . /t \ 1/4 



Sp(r,t) 



- -A - 

2y/2 n V t 



r(l/4)-2v^F z + 0(z 2 ) 



(44) 



whilst in the latter, a steepest descents analysis yields 

Bz 2 ^ 



5p(v,t) 



7r r 



exp[-(3/4)^ 4/3 ] . 



(45) 
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Both of the above results are in complete agreement with the scaling functions found by 
Brummelhuis and Hilhorst Q], from an exact lattice calculation for VMD on an infinite strip. 
This gives us strong confidence in the physical integrity of our continuum theory of VMD. 

For completeness we briefly describe the form of the mean OP density for arbitrary 
dimension d G [0, 2]. The scaling variable in this case is generalized to 



2(7) 



1-7 



(l-7)/2 



rA 



(46) 



(t/to) (1 ^ )/2 

Referring to Eq.((|l]), inverse FT and subsequent analysis yields the following results. For 
2(7) C 1 we find 

b r(i- 7 ) 



5p{r,t) 



(47r)^r( 7 + (l-7) 2 ) 



, f n \ 7(1-7) r 
A 2 7 to\ r 



1 + 0(z 2{1 -i>) 



(47) 



whilst a steepest descents analysis for z(j) 3> 1 reveals 



Sp(r, t) ~ exp 



1+7) 



2(7) 



2/(1+7) 



(48) 



This completes our study of the simplest VMD scenario; namely the effective diffusion 
of a tagged particle by a wandering vacancy. We have been able to give exact results for 
all dimensions. In the case of integer dimensions, our results are found to be in complete 
agreement with previous exact lattice studies. 



VI. EXTENSION TO A SLAVED FLUX-LINE 

In this penultimate section we shall consider a more complicated VMD scenario, both 
as an illustration of the utility of the coarse-grained approach, and also as a physical model 
for a VMD mechanism of flux lattice melting. 

We consider an anisotropic lattice consisting of well separated planes (in particular, one 
of the high-T c cuprates fll0|). Within each plane there exists a low density of wandering 
vacancies. We now imagine a flux line directed perpendicular to the planes, and strongly 
pinned by certain lattice impurities [TTJ] . If the binding energy is strong, thermal wandering 
of the line will be completely suppressed. However, a much weaker form of line wandering 
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may be driven by VMD of the pinning sites themselves (due to exchange with the low density 
planar vacancies, or vacancy aggregates). Since each plane has its own stock of vacancies, 
(which we disallow from hopping from plane to plane), the interactions between a given 
vacancy and the line are recurrent and we can expect slow and steady smearing of the mean 
density of the flux line. The physical existence and/or relevance of this mechanism deserves 
more detailed investigation (for instance, there are several competing pinning mechanisms 
within the material, one of which is actually due to the oxygen vacancies themselves IfLlH ). 

We shall describe this system by generalizing the continuum theory of VMD outlined in 
section II. First, we take for simplicity one vacancy within each plane. In the continuum, 
this is simply described by attaching a longitudinal coordinate to the vacancy position R. 
The equation of motion for the vacancies is then given by 



where the noise has zero mean and covariance (£ a (z, t)^(z', £')) = D5 atl 35(z — z')5(t — £'). 
The OP cf) now describes the probability density of the flux line. It is a function of a planar 
coordinate r, a longitudinal coordinate z, and time t. For a given z, the evolution of the 



- namely, an elastic interaction (which stems from the Josephson coupling between Cu-0 
planes ||12|| ). Thus the equation of motion for the OP is 



where v is an effective longitudinal elasticity. 

As an initial condition, the simplest choice is to take the flux line to be straight, and 
located at the origin: <j>(r, z, 0) = A5(r). Also, we could (artificially) start all the vacancies 
at the origin: R(z, 0) = 0. 

This model may now be analysed in precisely the same way as our original VMD model; 
namely through an infinite order perturbation expansion in powers of A. One must use an 
additional longitudinal FT to diagonalize the elastic coupling. The appearance of multi- 
ple longitudinal Green functions makes the explicit evaluation of the functions (xn) more 




(49) 



OP is given by Eq.([|), in two dimensions. The simplest longitudinal coupling shall be taken 



a, 




(50) 
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challenging than before. However, analytic progress seems possible and is currently being 
pursued. It is certainly of interest to calculate the //-dependence of the mean OP evolution, 
to see how effectively the flux line elasticity combats the driving forces of VMD. Preliminary 
results indicate that i) for v — ► oo the line fluctuations become gaussian and coincide exactly 
with the fluctuations of a single point within MFT (as described in section III), and ii) for 
v — > the line fluctuations have a singular dependence on the elasticity (and thus differ 
from the fluctuations of independent planar tags as described in sections IV and V). 

VII. CONCLUSIONS 

In this paper we have constructed and solved a continuum theory of vacancy-mediated 
diffusion. Our main intention has been to test the theory against exact results known 
from lattice studies [|]||. In particular we have thoroughly examined the evolution of the 
mean OP density in the case of a single vacancy smearing an originally sharply-peaked OP 
fluctuation (which corresponds to the lattice scenario of following the motion of a tagged 
particle due to vacancy exchange). At the level of mean field theory, we found the (almost 
obvious) length-time scaling for the OP density, but no sign of the non-gaussian fluctuations, 
which were the most interesting results obtained from the lattice studies. Therefore we 
pursued a more systematic treatment, based on an infinite order perturbation expansion. 
We presented an exact analysis of our theory in all dimensions, and complete agreement has 
been found with the dynamical scaling results obtained from the lattice. In particular, we 
find that for d > 2, the OP density is smeared over a limited range and then freezes. The 
envelope is a simple exponential. In d = 2, the smearing is slow, but continues indefinitely. 
The envelope is described by the modified Bessel function Kq. Finally, for d < 2, a more 
challenging calculation revealed that the envelope of OP smearing is described by a stretched 
exponential. Our results have the advantage of being valid for arbitrary dimension, thus 
revealing more clearly their analytic structure. In the final section we proposed an extension 
of simple VMD to the diffusion of a pinned flux line, slaved to planar vacancy exchange. 
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Analysis of this more challenging problem is in progress. 

There are, needless to say, many other extensions to the current work. We mention the 
more obvious here. First, it would be interesting to calculate higher order OP correlations. 
Throughout the present work, we have concentrated on the mean OP density. However, there 
is much non-trivial information hidden in the simplest spatio-temporal correlation functions. 
Second, one can apply the continuum theory to more complicated (single vacancy) scenarios; 
mainly by adjusting the boundary and initial conditions. For instance, a step function initial 
condition would correspond to the vacancy-mediated roughening of an initially straight 
domain wall. Finally, one could investigate simple mechanisms whereby the random walk of 
the vacancy itself is weakly coupled to the OP distribution, which is a physically relevant 
perturbation. 

The author would like to thank Z. Toroczkai and R. Zia for bringing this problem to his 
attention. The author also thanks E. Lundell for a critical reading of the manuscript. The 
author gratefully acknowledges financial support from the Division of Materials Research of 
the National Science Foundation. 
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APPENDIX A: 



In this appendix we explicitly evaluate the average (G(k, t)G(k 1 , t)*), where G is defined 
in Eq.(f^). Aside from momentum prefactors, we need to evaluate 



(Al) 



\ ' 

We have generalized the momentum in the exponent to a time- dependent form k(£) for a 
reason soon to become clear. At the end of the averaging procedure we shall reset k = k — kx 
as required. 

The average given above is most easily evaluated by generating the noise prefactors via 
functional differentiation with respect to n(t). Thus, 



d 2 



5K a (t)5^(t) 



exp 



exp 



i / dt' £(*') ■ K(t') 



D/2) j dt' K{t'f 



D 



Dt Q K a {t)K P {t) 



exp 



{D/2) / dt' n(t 



/\2 



(A2) 



where to is the implicit scale of the noise correlations. We need only retain the first term, 
given the smallness of to- Using this result with Eq.(|TTJ) yields Eq.(|P2"D in the main text. 

APPENDIX B: 

In this appendix we outline the evaluation of Q n as defined in Eq . (^2|) . The terms in 

the exponent of Q n are of the form R(t m ) • (k m _i — k m ) = / ds m £(s m ) • (k m _i — k m ). We 

o 

replace the momenta appearing in the exponent with generalized momenta n m (s m ). After 
the averaging procedure we set n m (s m ) = k m _i — k m . This allows us to generate the noise 
prefactors from functional differentiation. Thus we have 



Qn({^mi t m , tt m , Pm}) 



n 

.m=l 



exp 



n L J 

' J ds ™ Us 

m=l n 



ml ' ^m\Sm 



;bi) 
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The average appearing above is easily performed over the multivariate gaussian noise distri- 
bution, and yields 

1 2 



[exp 



m=l 



D n 

vx y> <i — 2 S 

m=l 



ds 



tm+l 



m 

.1=1 



(B2) 



A given double functional derivative of the exponent gives: 



8K%"(t m )8Kt(t m ) 



exp{-(D/2) J ds [K(s) + K m (s)} 2 } 



tm+l 



D 



x exp{-(£)/2) J ds [K(s) + /t m (s)] 2 } . (B3) 



tm+l 



We use this last result to perform the n double functional derivatives in Eq . (pl| ) . We note 



that using our final replacement for {ft m } yields — k — k m . Thus we reproduce Eq . (P3| 
as given in the main text. 



APPENDIX C: 



In this appendix we give a brief analysis of the integral appearing in Eq. ( f4T|) , specializing 
to d — 1. In terms of the density difference, this has the form 

1 



(CI) 



T 1 / 2 +b{k,t)_ 

where for convenience we have set b = (&;/A) 2 (t/to) 1//2 - First, we change variables from r 
to —x (which runs along the negative real axis) by integrating across the branch cut. This 
yields 

B °r dx e~ b2x 
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(C2) 



<x (1 + x) 

Next we change variables to y — \fx. Rewriting the exponential using a Hubbard- 
Stratonovich transformation gives the double integral 



dy 



5p(k,t) = — , 

n J (1 + 



ds 



-s 2 +2ibsy 



(C3) 
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The integral over y now resembles a Fourier transform, and results in a simple exponential, 
so that 

5p(k,t) = B [ ^e- 2 " 26 !-! . (C4) 

J V7T 

— 00 

The inverse Fourier transform from k to r is now easily performed resulting in Eq. (f4H|) , as 
shown in the main text. 

For general d G [0,2], the analysis is more difficult. Progress is made by exponenti- 
ating the Lorentzian form in Eq.(|4*ID using an auxiliary integral, and then performing the 
asymptotic expansions on the resulting double integral. 
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